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Abstract: I argue that the incoherent, zero-frequency limit of the universal crossover 
function in the temperature-dependent conductivity at the superconductor-insulator transi- 
tion in disordered systems may be understood as an analytic function of dimensionality of 
the system d, with a simple pole dX d = 1. Combining the exact result for the crossover 
function in d = 1 with the recursion relations in = 1 + e, the leading term in the Laurent 
series in the small parameter e for this quantity is computed for the systems of disordered 
bosons with short-range and with Coulomb interactions. The universal, low temperature, 
dc critical conductivity for the dirty boson system with Coulomb interactions in ci = 2 is 
estimated to be ~ 0.69(2e)^//i, in good agreement with many experiments on thin films. The 
next order correction is likely to somewhat increase the result, possibly bringing it closer to 
the self-dual value. 



Different quasi two-dimensional (2D) electronic systems, like tliin films Joseplison 
junction arrays or underdoped high-Tc cuprates 0, appear to have a continuous zero- 
temperature phase transition between the superconducting and the insulating state, as some 
parameter of the system is varied. The loss of phase coherence in the ground state is believed 
to be due to Anderson localization of the Cooper pairs 0, which through the quantum 
uncertainty relation competes with the phase ordering. In reality, this purely quantum 
(T = 0) phenomenon is unavoidably observed at low but finite temperatures, typically as a 
crossover from increasing to decreasing dc conductivity with temperature. Particularly in- 
teresting is the behavior of the conductivity near the critical value of the tuning parameter, 
which remains finite down to the lowest temperatures, and usually very close to (2e)^//i, 
the quantum unit of conductance for electron pairs. This near-universality of the metallic 
transport right at the superconductor-insulator (SI) transitions in 2D has been a subject of 
numerous theoretical and experimental investigations within the last decade. Scaling argu- 
ments imply that at T = and right at the quantum critical point the dc conductivity 
must indeed be universal, very much like the critical exponents. This insight has received 
a substantial amount of theoretical support in form of concrete numerical |]^, [|], ^ and 
analytical [0] , [|1^ , [|lT] evaluations of the T = critical conductivity for various universality 



classes of the SI transitions. The experimental evidence for the expected universality of the 
charge transport however, is still somewhat less convincing. In particular, measurements 
typically yield values of the critical conductivity two to three times larger than the calcu- 
lations Most of the theoretical studies [0, 0, 0, [|T^, [|ri|] of the critical conductivity 
however, are in significant conceptual discord with the experiment, as recently emphasized 
by Damle and Sachdev They argued that the uj = 0, T —>■ dc conductivity that is 
typically measured and the T = 0, u; — > conductivity that is usually calculated, originate in 
entirely different dissipation mechanisms, and while both should be universal, they have no 
reason to be equal. As an illustration, they obtained the conductivity in the hydrodynamic, 
incoherent transport regime huo/kBT — > at the simpler superfluid-Mott insulator (SF-MI) 
transition brought by an external periodic potential, near its upper critical dimension of 
c? = 3, and indeed found a different, and larger, value of the critical conductivity. 

In real disorder systems the insulating phase is presumably the result of Anderson local- 
ization and should therefore be a compressible Bose-glass (BG) 0, not the Mott insulator. 
It has been argued that, with disorder present, the limits ^ and T ^ may commute 
0, and the imaginary time Monte Carlo calculations indeed found little or no dependence of 
the critical conductivity at the SF-BG transition on the ratio oJn/T @, where uj^ is the Mat- 
subara frequency. Including disorder makes the problem difficult to study analytically, since 
then the transition seems to lack the upper critical dimension, or any other obvious limit 
in which at the criticality the system would be weakly coupled. Numerical techniques [0, 
P], on the other hand, determine the conductivity by extrapolating from the imaginary 
Matsubara frequencies, uJn = n27rkBT, and thus by its nature are not very sensitive to any 
structure it might have in the hydrodynamic regime. The understanding of the role of disor- 
der and the calculation of the experimentally measured low temperature dc conductivity at 
the SI critical point in 2D systems thus presents itself as a fundamental and unsolved prob- 
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lem. The purpose of this Letter is to propose its solution in form of a controlled expansion 
around the solvable case of the SF-BG transition in one dimension (ID). 

I will consider the system of disordered interacting bosons of charge e=K = 2e defined, 
for example, by the Bose-Hubbard, or Josephson junction array Hamiltonian with a random 
chemical potential [Q, |Tl|]. This model is known to possess a continuous SF-BG transition at 



T = 0, and should be appropriate for description of the SI transition observed in Josephson 
junction arrays, or ^He in random media. If the Coulomb interaction between bosons is 
added, it should represent the correct universality class for the transition in homogeneous 
and granular films. Although there is no upper critical dimension for the SF-BG transition 



in the usual sense [0, 0, [T^ the fact that the superfiuid phase in d = 1 and at T = 
exhibits only a power-law long-range order implies that d = 1 represents the lower critical 
dimension. Recently, this observation has been used by the author to formulate a controlled 
expansion of the universal quantities at the SF-BG transition in powers of a small parameter 
e = d — 1 In particular, and as argued below, in the limit T ^ the dc resistivity at 



the critical point is proportional to a certain power of temperature and to the value of the 
disorder parameter at the SF-BG critical point, which becomes infinitesimally small (~ e) 
as dimensionality of the system is reduced to d = 1. This suggests that the universal part 
of the low temperature dc conductivity at the critical point may be expressed as an analytic 
function of e, with a simple pole at e = 0. I compute the first term in the Laurent series for 
the real part of the critical dc conductivity around e = 0: 

.;,..o,r)^(^)^[-5;,iH-0(i)l|. (1) 

where c is a microscopic constant with units [lengthy /time, z = {d,l} is the dynamical 
exponent and x = {1,2} for the short-range and the Coulomb interaction between bosons, 
respectively. For d = 2 this leads to an estimate cxc ~ 0.69el/h for the Coulomb universality 
class, in agreement with many experiments on thin films [Q. A systematic perturbative 
procedure for a further calculation of the critical conductivity and of the critical exponents 
is outhned, and the sign of the next order correction in the Eq. (1) is discussed. I speculate 
that at the SF-BG critical point in 2D the real part of the conductivity is a continuously 
decreasing function of the ratio fiu) /ksT, quahtatively similar to the result m. d = 1 + e. 

In general, right at the critical point in a dimensional system the real part of conduc- 
tivity for the frequencies and the temperatures much smaller then some microscopic cutoff 
energy can be written as 

where Fd{x) is a universal crossover function. While I focus here on the SI transition, 
the above equation applies equally well at the critical points in quantum Hall systems ||T3[ 
and doped semiconductors |jT6[, for example. At d = 2 the power of the temperature- 



dependent dephasing length in the last equation vanishes, and the conductivity in units 
of el/h becomes a universal function of the ratio of quantum and thermal energy scales. 
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It then becomes evident that, in principle, at the critical point the conductivity in d = 2 
assumes two completely different values depending on the order of limits uj —>■ and T 0, 
that correspond to the values of the crossover function either at zero or at infinity, and 
measure completely incoherent or completely coherent transport, respectively. What is more 
surprising, although it arises as a natural consequence of the scaling law, is that even in the 
limit huo/ksT — ^ of incoherent, collision dominated transport, conductivity at the critical 
point in c? = 2 is still given by a universal number |]12| . To understand the dependence of 
the conductivity in this limit on dimensionality of the system, consider the effective low- 
energy action for the disordered superfluid that describes the SF-BG transition in the dirty 
Bose-Hubbard model in d = 1 fl^ , [|ll|] , H : 



n J Jo 

-DY, dTdT'cos2{e,{x,T) -ej{x,T')), (3) 

where (3 = h/ksT^ K is inversely proportional to the superfluid density, c is the velocity 
of low-energy phononic excitations, index i numerates the replicas introduced to average 
over disorder, and D is related to the width of the distribution of the random potential. 
The usual limit — > and a short-distance cutoff A^^ in the Eq. (3) are assumed. This 
effective action arises as the one dimensional realization of the density (dual) representation 



of the dirty Bose-Hubbard model at low energies [|TT], or as the bosonic representation of 
the disordered Luttinger liquid fl^. Invariance under a change of the cutoff implies that the 
conductivity in o? = 1 can be written in the scaling form 

.(.,r) = |L;(^-(6).cW.DW.^)|, (4) 

where K{b), c{b) and D{b) are the renormalized couplings at the new cutoff bA^^, and 
b = hcA/ksT. At low temperatures the conductivity is determined by the infrared stable 
fixed point of the scaling transformation. The result of the renormalization in the theory 



(3) for weak disorder is well known [0, [14|: under the change of cutoff the combination 
K, = Kc^, that is proportional to the compressibility, stays constant, K always increases, and 
small D is relevant for rj = Kc > 1/3, and irrelevant otherwise. There exists a separatrix in 
the 1] — D plane which ends in the SF-BG critical point at = 1/3 and D = 0. At u; = 



the scaling function in the Eq. (4) at weak disorder should behave as f ^ 1/ D{b) [T^, [|T 
and therefore right at the separatrix slowly (logarithmically, ~ (ln(6))^, [|I^]) diverges as 
& — cxo in (i = 1. Apart from the logarithmic correction that derives from disorder being 
dangerously irrelevant at the SF-BG criticality in ID, right at the transition the Eq. (4) is 
just a special case of the general scaling form (2), for d = z = 1. Since in = 1 + e disorder 



at the transition scales towards a finite fixed point value D(b) ^ ~ e, when b ^ oo [|T4 
the comparison of the two scaling laws in Eqs. (2) and (4) leads to the identification 

Fi+,(0) = /(1/3, 1, D% 0) = + 0(1), (5) 
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which expresses the central idea of this work. The leading term in the Laurent series for 
Ffi{0) is completely determined by the scaling function as in ID and by the infinitesimal 
value of disorder at the fixed point of the scaling transformation in d = 1 + e. 

In the remaining of the paper the complete function / near the SF-BG critical point 
in d = 1 is obtained, a new field-theoretic version of the recursion relations in d = 1 + e 
requisite for determination of the fixed point value of disorder is derived, and finally, the 
residuum at the pole at e = in the Eq. (5) for both short-range and Coulomb interactions 
between bosons is computed. 

The standard linear response formalism yields the conductivity in d = 1 

n,m=l 

Gnmi.^) is the temperature dependent, retarded, g = Green's function defined as 

= JdxJ^^ dre'^'-^rMx, r)^„(0, 0)),^„^^+,5, (7) 

where Tr is the standard time-ordering operator. The thermal Green's function in the Eq. 
(7) may be evaluated perturbatively in disorder and then analytically continued to real 



frequencies. Similar calculation has been performed before by Luther and Peschel [T^ for 
rj > 1, which corresponds to weak coupling in the equivalent ID fermionic system. The 
SF-BG transition at weak disorder is at 77 ~ 1/3, so I derive here a slightly improved version 
of their results which can be analytically continued into the transition region. Introduce 
the self-energy as G*j(a;„) = 6ij/{2{K/TT)u'^ + S*(ci;„)). To the lowest order in D it may be 
written as 

= 8D 1'^ dril - e^-"-)(T,e*29,(x,r)g-i29,(x,o)^^ _^ 0{D^), (8) 
Jo 

where the average is performed over the quadratic part of the action (3). The self-energy 
is thus itself a Green's function, which enables one to perform the analytic continuation to 
real frequencies by first rotating the integrand in (8) to real time hy r ^ it to find the real 
time time-ordered propagator, and from it finally to determine the retarded one by using 



the standard relation between them |20]. Performing the Fourier transform at the resulting 



expression then gives the retarded self-energy 



L (cj) = — - — ( )') sm(— je" / dt (9) 



cA hcK 2r] Jo (sinh(7rt))^ 

where C ~ 0.577 is the Euler's constant, and I assumed that hcA/kBT » 1, i.e. the 
continuum (low-temperature) limit. Appearance of a particular numerical constant is a 
consequence of the assumption that the dispersion is u = ck for all momenta < A; < A, 
and is a non-universal, short-distance feature. This nevertheless, does not compromise the 
universality of the conductivity at the transition, since any non-universal constant like C may 
at the end be absorbed into the definition of the running, dimensionless disorder coupling. 
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as will be done shortly. The remaining integral in (9) is convergent only for r] > 1, but once 
evaluated there exactly, may be defined via analytic continuation in the transition region 
1] ^ 1/3. Performing the integral, in the vicinity of 77 = 1/3 one obtains the conductivity in 
d = 1 to be 



iujc el 



""^""'^^ 7]{T)u^ + i2{kBT/hYW{T)g{huj/kBT) h' 

where W{T) = (Tc^D/c^A^e^'"){kBT/hcA)'n^^ is the dimensionless disorder variable, r]iT) = 
T] + W{T)n-^/^tan{n/2ri), and g{x) = (1 + {x / rr)^) tanh{x / 2) . Note that the result indeed 
may be cast into the scaling form as claimed in the Eq. (4). The self-energy acquired an 
imaginary part, which in the dc limit hu/kBT becomes proportional to temperature 
and to the temperature dependent disorder variable W{T). The significance of the point 
7] = 1/3 now becomes apparent: for 77 < 1/3 disorder variable W(T) scales towards zero with 
decreasing temperature, and the lowest order result in the Eq. (10) becomes asymptotically 
exact. At a finite frequency and at T = the real part of conductivity in d = 1 then becomes 
cr'(co') = (7rc/?7(0))5(ti;) and the system is an ideal conductor, li r] > 1/3 the perturbation 
theory breaks down, which indicates the entrance into the insulating phase. Notice that as 
rj — > l/3~ the coefficient in front of W{T) in the expansion for rj{T) becomes divergent, as it 
has a simple pole at r/ = 1/3. This is reminiscent of the dimensional regularization frequently 



employed in the studies of thermal critical phenomena [^, and indicates that the theory 
defined by the action in Eq. (3) becomes just renormalizable at 77 = 1/3. 

In the continuum limit, the effect of change of temperature on the low- frequency con- 
ductivity in d = 1 may be expressed entirely through the effective values of the coupling 
constants rjiT), W(T), and, if we had retained the momentum dependence of the propa- 
gator, the compressibility k(T). This, of course, is just the statement of renormalizability 
of the theory, with the temperature used as an infrared regulator. Taking into account the 
non-zero canonical dimensions of the coupling constants away from d = 1 fl^, the effective 
couplings close to r/ = 1/3 satisfy the differential equations: 

V{T) = z-\d - IMT) - :^W{T) + 0{W'{T)), (11) 

W{T) = {- -3)W{T) + 0{W\T)), (12) 
V 

k{T) = z^\d- z)k{T), (13) 

where x = dx/dhi{kBT/hcA). The d-dependent terms in Eqs. (11) and (13) may be inferred 
from the scaling of the superfluid density and the compressibility near the critical point P|, 
as discussed at length elsewhere |T^. Note that in the Eq. (13), unlike in the Eq. (11), 
there are no terms proportional to W{T). This is a consequence of the exact symmetry 
of the interaction term in the action (2) under 6'j(x,r) 6i{x,r) + h{x), for arbitrary 
function h{x), and fixes the value of dynamical exponent to z = d for the system with 
short-range interactions IQ. Linearization of the fiow close to the fixed point of Eqs. (11) 



and (12) gives the correlation length exponent v = (l/v3e) + 0(1), in agreement with 
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the result of the momentum-shell renormalization group [|T^]. The fixed point is located at 
W*{T) = TT^/^e/e + 0(e2) and r]*{T) = 1/3 + 0(e). Using the Eqs. (10) and (5) one then 
obtains the main result announced in the Eq. (1), for the short-range interactions between 
bosons. 

To make a comparison with experiments on thin films or on high-Tc cuprates 
one needs to include the long-range Coulomb repulsion between the electron pairs. A way 
to do this was proposed previously by the author in the ref. 14, where the long-range 
interactions was defined as V{r) = / d'^qexp{iq- r)/g'^~^, so that it coincides with the 
Coulomb interaction for d > 1, and with the short-range interaction precisely at d = 1. The 
calculation of conductivity in ci = 1 then remains the same, the only change now being the 
equation for the temperature dependent charge e^(T) instead for the compressibility in the 
Eq. (13), with {z — d) —>■ {z — 1), and consequently, {d — 1) —>■ {d — l)/2 in the Eq. (11) 



|T^. It then follows that with Coulomb interactions present z = 1, u = y2/3e + 0(1), and 

the fixed point value W*{T) = 7r^/^e/12 -|- O(e^). This yields the second result in the Eq. 
(1), for the Coulomb universality class. To the lowest order, the critical dc conductivity is 
larger if the Coulomb interaction is present. This may have been intuitively expected: a 
longer-range interaction suppresses the phase order more efficiently, so it takes less disorder 
to finally turn the system into an insulator. 

Although there is no very good agreement on the value of critical conductivity between 
the different experiments, most of the measurements [|l| on thin films are very close to cxc ~ 
l(2e)^//i, in a quite reasonable agreement with my lowest order estimate for the Coulomb 
universality class. On physical grounds, one may also expect the next order term in the 
Eq. (11) to be again negative [|14|, and that there is no second order correction in D in 



the imaginary part of the self-energy [IS]. The first expectation is based on the fact that 



disorder inhibits superfiuidity, and therefore should effectively increase the exponent (r/) in 
the power-law for the superfluid correlator at large distances. Also, in ID the coupling 
D measures the 2k p backscattering over the random potential in the equivalent fermionic 
problem, so corrections should renormalize only the forward scattering amplitude, which 
in ID is unrelated to conductivity. In that case the next order 0(1) correction in the 
Eq. (1) would be positive, presumably bringing the result closer to the experimental one. 
It is also interesting to note that for both universality classes the critical conductivity in 
the hydrodynamic regime obtained here turns out to be larger than the one in the coherent, 
T = limit ip, |11|. This also appears to be true for the SF-MI transition in a commensurate 



periodic potential |jT2[|. Since d = 1 represents the lower critical dimension for the SF-MI 
transition as well, a similar calculation to the present one could presumably be performed 
in that case, except that it would require a calculation of the second order in strength of the 
periodic potential. It would be interesting to compare the result obtained this way with the 
calculation near the upper critical dimension of the ref. 12. 

Even though the Eq. (10) has been derived here with the purpose of obtaining the 
universal conductivity in d = 1 + e, it is directly applicable to transport near the possible 
transition in two edges of the incompressible quantum liquid in a gated Hall bar [H]. Also, 
by continuity, the Eq. (10) implies that the crossover function Fdi^x) is a continuously 
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decreasing function of its argument for d = 1 + e, with a maximum ~ 1/e at a; = 0, and 
vanishing as ~ e/x'^^~'^^/^ for large x. Ahhough not ruled out, it does seem unlikely that in 
d = 2 this dependence on x = Tioj/kBT should completely disappear. In fact, the difference 
in estimated critical conductivities in completely coherent and incoherent regimes suggests 
that the situation in (i = 2 is most likely qualitatively similar to c? = 1 + e: the real part of 
conductivity should continuously decrease as a function of hu/kBT, interpolating between 
the two finite, u = and T = 0, limits. 

Finally, the field-theoretic formulation of the renormalization group transformation de- 
rived here has the advantage over the usual momentum-shell calculation [jl^, fl^ in that it 
facilitates a more systematic higher-order calculation. The observation that rj in d = 1 plays 
a role similar to dimensionality in the classical critical phenomena suggests a procedure anal- 
ogous to the standard dimensional regularization for the ci-independent part of the recursion 
relations. Adding the effect of dimensionality when ci > 1 as described in ref. 14 and as done 
in Eqs. (11) - (13) would then yield the higher order corrections for the exponent u and 
the critical dc conductivity. It would be very interesting to compare the results of such an 
analytical calculation with the experiments and the numerical simulations, as it could lead 
to a more definite understanding of the SF-BG quantum critical behavior. 

The author is grateful to NSERC of Canada and the Izaak Walton Killam foundation 
for the financial support. 
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